NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 1908 


GENERAL ALGEBRAIC METHOD APPLIED TO CONTROL 

ANALYSIS OF COMPLEX ENGINE TYPES 

By Aaron S. Boksenbom and Richard Hood 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 
July 1949 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 1908 

GENERAL ALGEBRAIC METHOD APPLIED TO CONTROL 
ANALYSIS OF COMPLEX ENGINE TYPES 
By Aaron S. Boksenbom and Richard Hood 


SUMMARY 

A general algebraic method of attack to the problem of con- 
trolling gas-turbine engines having any number of independent 
variables vas utilized employing operational functions to describe 
the assumed linear characteristics for the engine, the control, 
and the other units in the system. Matrices were used to describe 
the various units of the system, to form a combined system show- 
ing all effects, and to form a single condensed matrix showing the 
principal effects. This method directly led to the conditions on 
the control system for noninteraction so that any setting dis- 
turbance would affect only its corresponding controlled variable. 
The response-action characteristics were expressed in terms of the 
control system and the engine characteristics. The ideal control- 
system characteristics were explicitly determined in terms of any 
desired response action. 


INTRODUCTION 

The current development of gas-turbine engines indicates a 
future trend toward a wide variety of engine types. As new 
engines are developed or built by combining basic components of 
existing engines, the control problem presented by each engine 
type will be different. Control systems satisfactory for one 
engine type could not be expected to operate for another engine 
type. Even engines of the same type but of different size or 
operated under radically different conditions may have different 
engine characteristics that require special control character- 
istics. The current fluid and indeterminate state of gas-turbine 
development indicates a need for a general control analysis that 
could be applied to any of the new engine types as they are 
developed. 

The control problem for engines differs from that for many 
plants, servomechanisms, or electric networks, because the control 
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of more than one variable is desired. At least those variables 
that may cause failure should be controlled. For gas-turbine 
engines, failure may be caused by excessive temperature, speed, or 
torque. Exceeding specified burner operating limits may cause 
burner blow-out. Thus, the problem of controlling more than one 
variable is basic for gas turbines. 

In control design, limitations may be placed on such char- 
acteristics as time of response, maximum undershoot or overshoot, 
dead-band, or steady-state error. This problem is the subject of 
much of the control literature to date and various methods of exact, 
approximate, and graphical means of solution have been developed 
(reference 1). 

Controlling more than one variable generally introduces an 
interaction among the controlled variables. Making a new setting 
for one variable may cause, during the transient state, changes 
in the other controlled variables. If these other variables are 
operating at or near a maximum point, this interaction may cause 
excessive values and possible damage to the engine. It would 
therefore be desirable that each new setting of a controlled 
variable affect only the variable being set, thus giving separate 
noninteracting control of all the variable s being controlled. 

A general algebraic approach to control problems of multi- 
variable engines developed at the NACA Lewis laboratory is pre- 
sented herein. This method leads to the conditions on the control 
characteristics to eliminate interaction between control variables 
and to an explicit solution for control characteristics to yield 
any desired response action. 


SYMBOLS 

The following symbols are used in this report: 
a,b,c,e constants 

C or C(p) control function to which engine- dependent- variable 

errors are applied 

C' or C'(p) control function to which engine- independent-variable 
errors are applied 


E or E(p) 
E* 


engine-characteristic function 

square matrix using first i columns of engine 
matrix 
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F(p) 


i 

L or L(p) 

N 

N 

n 

P 

R or E(p) 


R' or R’(p) 

S or S(p) 

T 

T 

V 




X 

x 

x 

X 


determinant of E* matrix 

cofactor of element of E* matrix 

general operational function 

number of engine- dependent variables being controlled 
instrument-characteristic function 
engine speed 

measured value of engine speed 
number of engine- independent variables 
operational symbol 

response function for controlled engine -dependent 
variables 

response function for controlled engine- independent 
variables 

servo- characteristic function 
engine temperature 

measured value of engine temperature 

uncontrolled transient disturbance in engine- 
independent variable 

engine fuel flow 

measured value of engine fuel flow 

signal to engine fuel- flow servo 

setting of engine-independent variable 

engine -independent variable 

measured value of engine- independent variable 

signal to engine- independent variable servo 
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Y 

y 

y 

£ 

7 

T 

Subscripts: 


setting of engine-dependent variable 
engine -dependent variable 
measured value of engine-dependent variable 
propeller-blade angle 
signal to propeller-blade-angle servo 
engine- interaction factor (y = 
engine time constant 



j,k,r,t,v,u indices 

s setting of specific engine variable 


ANALYSIS 

The system considered herein is essentially that shown in 
figure 1. This block diagram describes a typical engine and con- 
trol system. The general system configuration chosen is the one 
employing negative feedback in which error in the controlled 
variable is fed back into the control. The engine is represented 
as a box with a characteristic function E the input of which is 
the engine- independent variable x and the output is the engine- 
dependent variable y. The characteristics of the ins trument 
measuring the. value of the controlled variable is denoted as L 
for which the input is the variable itself and the output is the 
measured value y. This error is the difference between the 
setting Y and the measured value of the variable y. The char- 
acteristics of the control are denoted as a function C for which 
the input is the error and the output is a signal x to a servo 
unit, which is the device that actually changes the engine- 
independent variable and completes the loop. The servo character- 
istic is denoted as the function S and may be considered as a 
part of the control or as a separate entity, as shown in figure 1. 

In addition to the intentional disturbance made on the con- 
trol setting, inadvertent transient disturbances may appear in 
any part of the system. One possible souree of such transient 
disturbances has been included as an uncontrolled change in 
engine- independent variable. The variable x is thus the sum 
of a controlled part and of an uncontrolled part V. 
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The characteristic functions used to describe the various 
units of the system are operational functions or Laplace transforms. 
The dynamic characteristics of gas-turbine engines, controls, and 
servos can be approximated by such operational functions at least 
for the small deviations of variables in .which the linear assumption 
is valid. These operational functions can be used without refer- 
ence to the initial conditions of operation before disturbances 
occur,, if the initial condition is a steady-state condition and the 
values of all variables are their deviations from the initial steady- 
state condition. In general, the unit can be described by a general 
operational or transfer function as follows: 

y = F(p) • x (1) 

The block diagram for this unit can be shown as a box containing 
the function F(p) for which the input is x and the output is y. 

One method of deriving the operational function for a unit 
that approximately follows a linear differential equation is as 
follows for the equation given by 

.. ' . 

y + ay + by » cx + ex (2) 

where the dots above symbols indicate derivatives with respect to 
time. This equation can be algebraically represented in terms of 
the differential^ operator or in terms of the Laplace transform 
variable (reference 2). If the differential operator p is sub- 
stituted fer — , equation 2 can be written 
dt 

p 2 y + apy + by = cpx + ex (3) 

If this equation is solved for y, an explicit expression for the 
output can then be obtained as a function of p operating on the 
input x as follows: 

y = - c -E- t 1 . x - p(p) .x • (4) 

p^ + ap + b 

For each unit of the system, an equation similar to equation (4) 
can be written in terms of the input, the transfer function, and 
the output. A set of simultaneous algebraic equations is then 
obtained that can be combined or manipulated in any desired manner. 
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The resulting transfer functions can then be interpreted from the 
domain of the variable p to the real time domain as time responses 
of a variable to some specified disturbance. 

For the system of figure 1, the following equations can be 
written: 


Engine characteristic: 

y = E • x 

Control characteristic: / 

x = C • (Y-y) 

Instrument characteristic: 

Servo and transient disturbance : 

x = S • x + V 




J 


(5) 


Equations (5) can be combined for the effect of the disturbances Y 
and V on the controlled variable y as follows: 


= _MC . y + S . v 

ESCL + 1 ESCL + 


( 6 ) 


The previous procedure is, in general, followed herein. Equa- 
tions, such as equations (5), are written for each unit in the 
system and are then combined, as in equation (6), to obtain the 
effects of the disturbances. 


General representation of engine, - If the engine has more 
than one independent variable, there are additional degrees of 
freedom in control. The engine- characteristic function must be 
expanded to show how each dependent variable of interest is 
affected by all the engine inputs as follows: 


y l = E U • *1 + *12 • *2 ♦ 
y 2 = ®21 ’ X L + e 22 * x 2 + 

• • • 

• • • 

y i = E il * X 1 + E i2 ’ x 2 + 


+ E ln * *n 
+ E 2n * % 

+ E in * ^ 


> ( 7 ) 

J 
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In equations (7) each dependent variable y is expressed as a 
linear Siam of effects due to each independent variable x. For n 
independent variables, there are n terms in each equation. There 
is an equation for each dependent variable of interest. For i 
dependent variables, where i % n, there are i equations. 

The number of degrees of freedom of the entire system cannot 
be greater than the number of degrees of freedom of the engine 
alone. For continuously acting controls (excluding such controls 
as limiters), the total number of engine variables being controlled 
should not be greater than the number of engine- independent vari- 
ables. Equations (7) can be written in the following index form: 

n 

y J ■ Y, ^ (8) 

k=l 

where 

j = 1,2, . . . i 

The functions Ej k represent engine characteristics whereby each 
input separately affects each output. 

Equations (7) or (8) can be visualized in the matrix form of 
figure 2. The engine functions are set up in a rectangular array 
where Ej k represents the element in the j^h row and the k tJl 

column. The x inputs enter the matrix in the columns and the 
y outputs are attached to the rows. Each input multiplies all the 
elements in its column and each output is the linear sum of the. 
resulting products in its row. Thus the effects shown by equa- 
tions (7) and (8) can be pictorially traced in figure 2. For later 
use, a square matrix obtained by using only the first i columns 
in figure 2 is denoted by E*. 

General representation of control. - The control system is 
generalized to this extent: (l) Inasmuch as the engine has n 

independent variables, there are n degrees of freedom for the 
entire system. If only a limited number i of engine- dependent 
variables are to be controlled, the difference n-i of the engine- 
independent variables can be controlled. (2) The control system 
employs negative feedback where only errors are applied to the 
control. (3) Each error is to affect every engine- independent 
variable. 
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The control equations are written where each control output x 
is expressed as a linear sum of effects due to the errors in the 
controlled variables as follows: 


x x = C n • (Y^-y^ + C 12 • (X 2 -y 2 ) + 

+ C ’l(i+1) ’ ^ X i+l" x i+l^ + * 
5n = C nl * ’(Vyi) + C n2 * < Y 2 + 
+ c 'n(i+l) * ( x i+l-^i+l) + • 


■'ll 


+ C’ 


(Xi- yi ) 


A 


In 


+ C. 


ni 


( Y i-Xi) 


>(9) 


+ C'nn • (Xn-xjJ 


There are a total of n controlled quantities, n errors, and 
n terms in each equation. There are also n engine- independent 
variables and therefore n equations. Those independent variables 
to be controlled are given subscript numbers from i+1 to n and 
those controls to which the corresponding errors are applied are 
denoted by C'. 


The control characteristics of equations (9) can be written 
in the following equation in index form: 


where 


i 



v=l 


n 



U=i+1 


<vv 


( 10 ) 


k = 1,2, . . . n 


Equations (9) and (10) can be visualized in the matrix form of 
figure 3 . The control functions are set up in a square arrange- 
ment where represents the element in the k^S row and 

v th column. The error inputs, Y-y and X-x, enter the matrix 
in the columns and the x outputs are attached to the rows. Each 
input multiplies all the elements in its column and each output is 
the linear sum of the resulting products in its row. Thus, the 
effects shown by equations (9) and (10) can be pictorially traced 
in figure 3. 

Representation for Instrument and servo functions. - The char- 
acteristics of the instrument measuring the value of the controlled 
variable has been denoted as L in equations (5), the input of 
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which is the variable itself and the output is the measured value. 
These characteristics can be written as 


where 


where 



I «w * y v 




v = 1,2, . . . i 

■ V • ^ 


(id 


\i = i+1 . 





The general assumption is made that each engine- independent 
variable would require a servo, which would vary the variable in 
accordance with its corresponding signal x. A source of transient 
disturbance is included that may vary any of the engine -independent 
variables outside of the controlled variation in these variables. 
Thus, the expressions for the engine- independent variables can be 
expressed by 

- Sfck • ik ♦ \ 

where 

i 

k = 1,2, . . . n (12) 

Each variable x is thus the sum of a controlled part and an uncon- 
trolled part. 

Analysis of complete system. - The determining characteristic 
equations for all the units in the system are repeated. 


\ 
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Engine: 


y j = 7 E Jk * x k 


( 8 ) 


where 


k=l 


3 = 1,2, . * . i 


Control: 


x 11 

5k ^kv • (Yy-y v ) + ^ ^*lqi * (^p -2 p) (10) 


where 


Instrument: 


where 


where 


V=1 n=i+l 

\ 

k = 1,2, • • • n 

y v = * yy 

v = 1,2, . . • i 

* V • \ 

(i = i+1 . . . n 




>(13) 


> (ID 


J 


Servo and transient: 

*k = Sjjt * 5k + v k (12) 

where 

k ~ 1,2, • • • n ^ 

These equations can he visualized in the matrix form of, 
figure 4. This figure represents, by block diagram, the general 
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system considered herein. This figure specifically is for an 
engine having five independent variables for which three dependent 
and two independent engine variables are controlled. The engine 
matrix is shown for the five inputs and the three outputs to be 
controlled. The control matrix is shown for the three y error 
Inputs, the two x error inputs, and the five outputs to the 
engine servos. The measuring devices, the servos, and the tran- 
sients are also included. The entire system is enclosed except 
for the outside disturbances, which can be imposed on the system. 
These are the five settings of controlled variables Y and X 
and the five transient disturbances that have been included. 


All the properties of the system can be obtained from this 
block-diagram representation. The effect of any outside disturb- 
ance or any internal disturbance can be easily traced on figure 4. 
For specific effects, the equations of each unit of the system, 
equations (13), are most productive of results. 

When equations (13) are combined to give the effects on the 
controlled variables of the outside disturbances the following 
equations are obtained: 


n i 


r J = C kv ^ Y v _I Vv y v^ + 


k=l v=l 
n n 


^ Stt cv <VV v - 


k=l p=i+l 


where 


E Jk V k 


k=l 


( 14 ) 


j — 1*2, . . . 


i 
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j. n 

x k ®kk ^kv(^v"^*vyyv^ + ^ ®kk^ , ku(^u“* J uu *u) + ^k 

V =1 U=i+1 (15) 


where 


k s 1,2, • • • n 


Equation (15) is valid for the effects on any (k = 1,2, . . . n) 
independent engine variable. For the effects on the controlled 
independent engine variables, equation (15) is used for k = i+1 
to n. A more condensed form of the matrix representation of the 
system, shown in figure 5, is indicated by equations (14) and (15). 
A single-system matrix is used in which the inputs are the errors 
in controlled variables and the outputs are the controlled vari- 
ables. The ESC matrix is a matrix in which the element in 

n 


row and v^* 1 column is 


I 

k=1 


E jk Skk C kv* 


Similarly, an ESC' matrix 


is 

is 


a matrix in which the element in the row and 


E 


jk Skk c 'k|i* An SC matrix is defined as Cj^y 


column 
and an 


k=l 


SC' matrix is defined as Sj^ C'j^. An additional matrix is 
included showing the effects of the transient disturbances. 


Noninteraction conditions . - Several types of noninteraction 
may be specified as follows: 


(1) A specific controlled variable, such as yj, is to be 

affected only by its own setting Yj and is not to be affected 
by any other setting, which means that in figure 5 the row 

of the system matrix contains only the diagonal element the 

other elements must be zero. 

(2) Any specific setting, such as Y t , is to affect only its 
corresponding controlled variable y t and not any of the other 
controlled variables. In figure 5 the t th column Of the system 
matrix contains only the diagonal element and the other elements 
must be zero'. 
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(3) Every controlled variable is to be affected by its corre- 
sponding setting only or every setting is to affect its correspond- 
ing controlled variable only. This complete noninteraction means 
that in figure 5 the system matrix is a diagonal matrix. In gen- 
eral, the effect of the setting on the system is of primary inter- 
est and a condition such as (1) is of secondary interest. The 
interaction effects defined in terms of the effects of the settings, 
such as (2) and (3), will thus be considered herein. 

For any setting such as Y t to affect only y t , it is suf- 
ficient from figure 5 that for all j from 1 to i vhen J £ t 


n 

E jk S kk C kt = 0 

k=l 

and for all p from i+1 to n 

Spp^pt “ ® 


(16) 


(17) 


or 


xt. 

Equation (17) states that the i+1 to n elements of t column 
of the C matrix (fig. 3) are to be zeroed. Along with this con- 
dition, equation (16) becomes 



E Jk s kk C kt = 0 


k=l 


(18) 


for all j from 1 to i where j £ t. Equation (18) expresses 
i-1 homogeneous equations in the i unknown C's and therefore 
specifies the ratios between any two elements in the t"“* column 
of the C matrix. This result becomes 


C Jt = 
®vv ®vt 


E* 

E* 


si i 

tv 


(19) 



14 


MCA TN 1908 


as shown in the appendix. The notation E* refers to the square 
engine matrix using only the first i columns. The notation 

indicates the cofactor of the Rj-j element of the E* matrix. 

The necessary and sufficient conditions for any setting to 

affect only y t are 

c pt = 0 (17a) 


for all u from i +1 to n and 


S , 1,1 C . 1 t = j^jyJ 

S w c vt tv | 


(19) 


as shown in the appendix. 

Now, for any setting X r to affect only Xj. it is sufficient 
from figure 5 that 

n 

E jk S kk c 'kr = 0 ( 20 )' 

k=l 


for all j from 1 to i and for all 4 from i +1 to n where 

S 44 C Vr = 0 (21) 

or 

I 


Equation (21) states that the i +1 to n elements of the r^ 
column of the C’ matrix contains only the diagonal term. With 
this condition, equation ( 20 ) becomes 

1 

^ (®jk ®kk ^ kr^ + ®jr ®rr ^*rr = ® (22) 


k=l 
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for all j from 1 to i. Equation (22) expresses i homo- 
geneous equations in the i+1 unknown (C* ) >s and therefore 
specifies the ratios between any two elements in the r^ col umn 
of the C* matrix. This result becomes 


t 


E* 


S jj 01 Jr _ kfl_ 

SC' i 

w vr 


kj 


®kr 


I 

k=l 


E* 


kv 


Ekr 


(23) 


for any j find v from 1 to i or 

i 

sc * r \**i 

= k=l 1 


hx 


S C' 

rr rr 




(24) 


as shown in the appendix. The necessary and sufficient condition 
for any setting Xj. to affect only Xj. are: 


C V - 0 

for all |ji from i+1 to n when ^r and 

i 


(21a) 


-£ 


I E*. 


S .1,1 C> Jr = k=l 


kj 


Ekr 


S™ C' _ 
rr rr 


IE* 


(24) 


as shown in the appendix. 

For complete noninteraction, where every setting is to affect 
its corresponding variable only, the conditions can be summarized 
as follows: 


Sit -° 


(17a) 


where 


p. = i+1 . . . n 
t = 1,2, . . . 


i 
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where 


where 


where 


C' 


nr 


= 0 


n = i+l . . . n (n^r) 

r = i+l . . . n 


S JJ C Jt = 
^vt 




(21a) 


(19) 


(24) 


J — 1,2, • • • i 
r = i+l .. . . n 

Equations (17a) and (21a) state that the i+l to n rows of the 
control matrix contain only diagonal elements. Equation (19) gives 
the required ratios between any two elements of any column of the 
C matrix. Equation (24) gives the required ratios in any column 
of the C' matrix of any element to the diagonal element. < 

Besponse equations. The use of the operational functions to 
characterize each unit in the system allows algebraic representa- 
tions of the relations among the variables of the system. With 
complete noninteraction specified, the system matrix of figure 5 
becomes a diagonal matrix and each circuit acts as follows: 

- n 

V E J* v k <25) 

k=l 


k=l 
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and 


*j. ~ Spr C * ) + 7p (26) 

When these equations are solved for the controlled variables as a 
function of the settings and the transient disturbances , the fol- 
lowing equations are obtained: 


E jk hk °kj 


7a = 


k=l 



E jk V k 


3 ~ i 


* t J*T 


k=l 


( E jk S kk °kj + 1 / ( E jk s kk c kJ L jj) + 1 


k=l 


k=l 


(27) 


s„ = 


S w C* 
rr rr 


^rr ^'rr ^rr + ^ 


s rr ^'rr *rr + ^ 


(28) 


If the operational response function of yj to Yj is defined as 
Ej j and the operational response function of Zj. to Xj. as 
E' rr , equations (27) and (28) become 



n 



l 33 - 1 ) E Jk v k 


k=l 


(29) 


Xy. - . Xr - (B*„ Lrr - 1) V r (30) 


where 
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E jk hk 


E<< = 


k=l 


JJ " i 



( E Jk ^ °kj L jj) + 1 


k=l 


(31) 


and 


B' 


S„ C’_. 
rr rr 


rr 


Spr V'rr ^ rr + ^ 


(32) 


The response of any yj to Yj insofar as the control matrix 
is concerned depends only on the elements in the column of 

the control matrix (fig. 3). Inasmuch as for complete noninteraction 
only the ratios of the elements of that column of the control matrix 
are specified, the response can be independently set by choosing any 
one element of that column. If the conditions for noninteraction are 
used, Rjj is a function of any C v j. 


JJ |E*| S„ Ljj + 


(33) 


Equations (32) and (33) show the response of any controlled variable 
to its setting as determined by the control functions. The response 
of a controlled-engine-independent variable does not depend on 
engine characteristics (equation (32)). The response of a controlled- 
engine- dependent variable depends on the square engine matrix E* 

(fig. 2). 

An important result is obtained when equations similar to (29) 
and (30) are written for the errors in controlled quantities as 
follows : 


T r y j - (1 - e jj ( Y J - 



k=l 


(34 \ 
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= " ®*rr * * r r) (^r ” ^ *r r ^r) (35) 

The errors in controlled quantities respond in the same manner 
to both setting and transient disturbances and the control con- 
figuration chosen can act both as a control to setting disturbances 
and a regulator to transient disturbances, at least with the tran- 
sient disturbances assumed. 

Ideal control functions. - The control-design problem is one 
of fitting a control system to a particular plant or engine in 
order to control certain variables in seme specified manner. The 
nature of the final system is characterized by the desired response 
of controlled quantities to outside disturbances. These responses 
may be exactly specified or given in general qualitative terms. 

The nature of the engine or plant to be controlled must be known 
in its static and dynamic states. Any method of explicitly solving 
exactly for the control functions for any given plant and corre- 
sponding desired response action would give an exact answer to this 
problem. 

The use of the operational functions to describe the engine, 
control, and response functions and the algebraic methods employed 
herein leads directly to an explicit solution for the ideal control 
functions. Equations (32) and (33), which give the response actions 
in terms of the control functions, can be solved for the control 
functions in terms of the response functions as follows: 



It is therefore theoretically possible to choose any desired 
response characteristics and to solve for the required control 
functions to achieve this response. The use of these equations as 
well as the noninteraction conditions previously given for control 
design have certain limitations that preclude the use of this 
method in its most general sense. These limitations arise from 
the inability of exactly incorporating in a control the engine 
characteristics that the control is to match. A simple first-order 
response can be called for in equation (36) in spite of the higher- 
order transfer functions of engine and servo. This type of 
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response is possible only if the control function can reduce the 
order of the transfer function by allowing factoring of the numer- 
ator and denominator of equation (33) into exactly the same terms. 

The control functions found from equations (36) and (37) for 
responses that are very stable would be expected to give satis- 
factory results for small changes in characteristics and in the 
engine characteristics the control function is to match. Special 
cases of engine characteristics may require that an unstabilizing 
term in the engine functions be factored out by the control. In 
this case, variations in engine or control characteristics would 
unstabilize the system. In such cases, equations (36) and (37) 
should be used with the proviso that the control does not decrease 
the order of the response function and the response functions must 
be chosen accordingly. 

The algebraic methods used in dealing with operational func- 
tions can be further generalized in terms of equivalent relations 
when the engine, the control, and the other unit characteristics 
are expressed in other forms. For instance, if the frequency- 
response characteristics are graphically known from tests or anal- 
ysis, all the results described herein can be graphically performed 
in terms of the frequency-response characteristics in place of the 
operational functions. 

Summary of basic theory and results. - For an engine having 
n engine- independent variables and for which i engine- dependent 
variables are to be set and controlled (i < n) and for which 
n - i engine- independent variables are to be set and controlled, 
engine, control system, instrument characteristics, servo, and 
transient disturbances are 

Engine: 

. n 

y J ■ < 8 > 


where 


j = 1,2, . . . 


i 
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Control : ± 


n 



2k =/ C kv (Vy v ) + > C V * (10) 


v=l 


where 


Instruments: 


H=i+1 


k = 1,2, • • • n 


*v ■ L w • 


where 


where 


V s V *4 

v = 1,2, t « • i 
|i = i+1 . . . n 

Servo and transient disturbance: 

x k = S kk * 2k + ^k 

k = 1,2, • • • n 


( 11 ) 


( 12 ) 


The necessary and sufficient conditions for any setting Y t 
to affect only y^. are 


c ut -° 


(17a) 


for all u from i+1 to n 


S .1.1 C Jt *tj 


*Vt ®*tv 


(19) 


The necessary and sufficient conditions for any setting Xp 
to affect only Xp are 
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C* =0 

\1T 


for all \i from 1+1 to n vhen \i£r 




rr 



( 21 ) 


(24) 


For complete noninteraction, equations (17) and (19) must be true 
for all values of t from 1 to i and equations (21) and (24) 
must be true for all values of r from i+1 to n. 


With complete noninteraction specified, the response of the 
controlled variables and their measured error signals to setting 
and transient disturbances are 


and 



*r - B 'rr • - < E ’rr T r 


= L JJ* 





(29) 


(30) 


(34) 


rr (^r^^rr (35) 


where 


* 
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= “T 


^~^ E Jk ^ C kj 
k=l 


M S W C V,j 


( E jk ^ c kj L JJ ) + 1 


H S w C vJ Ljj + ^^1 


k=l 


(31) and (33) 


and 


E* 


S__ C* 
rr rr 


rr 


^rr C'rr **rr + ^ 


(32) 


With complete noninteraction specified, the explicit expression 
for the ideal control functions in terms of the desired response 
functions are 


E 


11 


“V ° T J = TT - *jj Ljj) 



(36) 


E' 


s rr C’ rr 


rr 


1 - E' 


rr 


krr 


(37) 


SPECIFIC EXAMPLE 

If the turhine-propeller engine is considered, which has the 
two independent variables, fuel flow and blade angle, the defining 
differential equations obtained by linearizing the general func- 
tional forms of engine torque as a function of speed and tempera- 
ture, propeller torque as a function of speed and blade angle,' and 
temperature as a function of speed and fuel flow are 

o 

N + TN = - a8 + bw f (38) 

T = cwj - eN (39) 

When these equations are placed in operational form and speed 
and temperature are solved for in terms of fuel flow and blade 
angle, 
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I 


N = 


1 + Tp 


P + 


1 + Tp 


(40) 


T = 


ae 


1 + Tp 


. p + ctoi. W£ 


(41) 


where 



' Case I 

i 

Controlling N and Wf. - For the case where n=2 and 1=1 
the engine matrix of figure 2 is obtained from equation (40) and 


E 


11 


-a 

1 + Tp 


E 


12 ~ 1 + Tp 


= P x 2 = v f 

y x = n v 

The general control system is 

P = C U (N S -N) + C* 12 (w fs -w f ) 

w f = C 21 (N S -N) + C* 22 (w fs -w f ) 

From equation (17a), the conditions for complete noninteraction 
give 


and from equation (24) for complete noninteraction. 


S 11 C *12 
S 22 C *22 


E*n E 12 


E* 


12 

11 


b 

a 
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where the determinant of a single element matrix written as 
|E»| = p]jJ E 1x allows the use of the equations used herein 

for this case. 

The control- system equations thus become 

g 

l . C u (N e -N) + -| jg C ' 22 (w fs -»f) 
if ’ °’22 

The responses of speed and fuel flow to their settings become 
from equation (33) 



E* 

s n c n 


E n s u c n 

E* S 11 C 11 L U + 

■*ii 

E n s n c n Hi + 1 


a S 11 C 11 

a S-q C-^ - (l+Tp) 

and from equation (32) 


, S 22 E '22 
22 S 22 C 22 h > 2 + 1 

\ 

The expressions for the ideal control functions then become from 
equation (36) 


B ll- | **11 1 R 11 1 _ _E 11 (l+Tp) 

'll 11 " 1-E]^! lax p| (1-E U Ln) E u (l-E u 1^7 a 


E’ 


22 


S 22 C ' 22 1 - B’22 L 22 


and from equation (37) 
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Case II 

Controlling N and T. - For the case vhere n=2 and 1=2, 
the engine matrix is obtained from equations (40) and (41) and 


E-n = - t — 2 

11 1 + Tp 

E 12 = 1 + T p 

•ni _ a© 

J c(i-r+rp) 

21 1 + Tp 

E 22 - 1 + Tp 

^ - P 

II 

K 3 

II 

tH 

II 


The general control system is 

\ 

P = C 11 (N S -K)+ C 12 (T s -T) 

v f = c 21 (n s -n) + c 22 (t 6 -t) 


The conditions for complete noninteraction give from equation (19) 


S 11 C 11 b 3*11 
S 22 C 21 | E *‘l2 



- 5: d-rtTp) 


S 11 C 12 
S 22 C 22 


E* 

E* 


Hi 

22 


E 12 _ b 
E 11 a 


The control- system equations thus become 


1 » c n<V B > + 1 J 2 c 22 (t 6 -t). 


11 


* * t(l-Arp) C 11<V H > * <WV T > 


The responses of speed and fuel flow to their settings become 
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, M s n c n s n c n 

11 ~ |®*| S 11 c ll Hi + l^lll ~ S 1X Cu. Igx - -IH7+ T P? 

61 


or 


Ell = 


S 22 C 21 


S 22 C 


21 Hi + b 


B 


**l S 22 C 22 


22 IE* I S 00 C_„ L 00 + IE* 


S 22 C 22 


22 22 22 |" 22 S 2 2 C 22 1^2 + c 


The expressions for the ideal control functions then became 


Sn C 


Hi 


IE* 


11 


- E 


ii (j-7+Tp) 


11 ■ 1 - En L n |E* | = 1 - Eu Lu a 


O r R 22 1^221 B 22 1 

22 22 1 - Egg Lgg |E*| (l " K 22 I^g) c 


SUMMAEY OF RESULTS 

The control of gas-turbine engines was limited to the linear 
assumptions for the engine, control system, and other units in the 
system. Operational functions were used for these units, the sys- 
tem was generalized to allow for complex engine types having any 
degree of freedom, and algebraic methods were employed in setting 
up the characteristics of the engine, the control, the servo, 
the measuring instruments. These algebraic methods were employed 
to combine the various separate pnits into a system showing how 
various effects are transmitted throughout the entire system. A 
condensed matrix form was used to illustrate the principal effects 
on controlled variables. 

By use of either the index or matrix form of the algebraic 
equations of the system, the following principal conditions were 
explicitly expressed: 

(a) The conditions that must be Imposed on the control func- 
tions so that any setting disturbance will affect only 
. the variable being set 
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(b) The conditions that must be imposed on the control func- 

tions for complete noninteraction so that every setting 
disturbance will affect only its corresponding variable 

(c) The response functions, which characterize the controlled- 

systern response action, in terms of the control- system 
functions 

(d) The ideal control-system characteristics in terms of any 

desired response action 

A specific example using the turbine-propeller engine having 
two degrees of freedom was given to illustrate the use of the 
previous results* 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 25, 1949. 
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APPENDIX - CALCULATION CONDITIONS 
Conditions on C Matrix 

Necessary condition. - The general equations for the response 
of the controlled variables are given by equations (Al) and (A2), 
(equations (14) and (15), respectively, in the text) as follows: 

n i 

y j ^ ^ C kv + 

k=l v=l 

n n n 

Y" «jk ^ ° v < w v *2^ Ejk Tk (A1) 

k=l \i= i+1 k=l 

where 

j = 1,2, . . . i 



x k “ / ®kk ^kv ^ y v“^w y v ) + ®kk 9 ]qi + ^k 



v=l 


U=i+1 


where 


(A2) 


k = 1,2, . . . n 

If the effects of any setting are separated, equations (Al) 

and (A2) become 

n i 

y j = /, /Y E jk S kk C kv <V L w y v^ + 
k=l v=l 
v/t 

n n 

IE 0, b. <w V ♦ 

k=l n=i+l 
n * 

®jk-Yk + / , E jk ®kk ®kt ( y t"**tt yt) 
k=l k=l 


. (A3) 
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i n 

Ckv (^y-^w ^y) + ^ ®kk ^ kn 

V=1 \i=±+l 

v^t 


v k + ^ c kt ( Y t ~ ^tt yt) 


(A4) 


If any setting Y^. is to affect only y^, then any yj is not 
affected by Y t ( J^t) and any controlled x k (k=i+l to n) is not 
affected by Y t . If equation (A3) is considered for any $£t, then 
the only term that can be a function of Y-t is the last one. 
Because of the arbitrary nature of the instrument function, 

Y t^ L ttYt 80 that 



E jk Skk C kt ■ 0 


(AS) 


k=l 


where 


J t t 


If equation (A4) is considered for k=i+l to n, the only term 
that can be a function of Y^ is the last one. Therefore, 


Skk c kt s 0 


c kt -o 


(A6) 

(A7) 


where 


k = i + 1 to n 

If equation (A7) is used, equation (A5) becomes 


i i 

^ N ® jk ®kk ^kt ® jt ®tk ^kk ^kt 

k=l ' k=l 
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Because j is any number (1,2, . . . i ), a similar expression 
for any other is 



Syy Cyfc - 


|^tv| y ^tk ^kk ^kt 

L 

k=l 
|E»| 


(A9) 


The ratio of equations (A8) and (A9) becomes 


S JJ C jt = l^tj 


Syy Cy-(. ^tV 


(A10) 


where 


j = 1,2, . . . i 
v = 1,2, . . . i 

Equations (A7) and (A10) thus follow necessarily from the noninter- 
action condition specified. 

Sufficient condition. - If equation (A10) holds, 

i 


r J ^VV C v t 



J=1 

where 

r = 1,2, . . . i 

From the theory of determinants, the expression 

i 


®rJ 


E* 


tj 


= 0 


J=1 
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when 


r ^ t 



when 


Therefore, 


r = t 


E rj S JJ c ,)t - 0 


when ' 

r^t 

If equation (A7) also holds, then equation (All) becomes 


n 

/ C Jt = 0 

j=l 

where 


(All) 


(AI2) 


r^t 

r = 1,2, . . . i 

If equations (A12) and (A7) hold, then any setting Yj. will affect 
only y t (fig. 5). Equations (A7) and (A10) thus are sufficient 
for the noninteraction condition specified. 
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Conditions on C' Matrix 

1 

Necessary condition. - When the general equations for the 
response of the controlled variables (equations (Al) and (A2)) are 
written separating the effects of any setting X r , these equations 
become 


n i 


n n 


r J Jk s kk c kv( E jk^ C 


k=l v=l 


k=l u=i+l 

nA* 



E jk v k + / E Jk s kk c W^~ L rr x r) (A13) 


k=l k=l 


JL 

■r- 



x k = y? ®kk ^kv (VWv) + ^ ®kk ^'k|i ) +• 

v=l ^=i+l 

p^r 


^k + ®kk ^'kr (^r"^ , rr x r) 


(A14) 


If any setting Xp is to affect only Xp, then any y^ is not 
affected by X r and any controlled x k (k=i+l to n) is not 
affected by Xj. for k/r . In equation (A13) the only term that 
can be a function of X r is the last one. Because of the arbi- 
trary nature of the instrument function, Xj. f L rr x r and 




jk Skk C 'kr 


= 0 


k=l 


(A15) 


/ 


where 


J = 1,2, . 


. . i 
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If equation (A14) is considered for any k(i+l to n) where k^r, 
the only term that can be a function of Xj. is the last one and 
therefore 

Skk C'kr = 0 (A16) 

and 


C ‘kr " 0 


(A17) 


where 


k/r 


If equation (A17) is used, equation (A15) becomes 
i 

"\ 

E jk ®kk ^'kr = ~ E jr ®rr ^'rr 



k=l 


where 


J = 1 . . . i 

The solution of this set of i equations in the i unknown 
Skk c 'kr e ives 


E* , E„ s C' 
pj pr rr rr 


q pi _ P = l 

C Jr 


E*| 


or 



(A18) 
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where 


J ~~ 1,2, • • • 1 

Equations (A17) and (A18) thus follow necessarily from the non- 
interaction condition specified. 

Sufficient condition. - If equation (A18) holds, 

i i 

S 

i 







S 3J C 'jr p-1 j=l 
s rr c ’rr |E*| 


= 0 


From the theory of determinants applied to the second term 

i 

C' 

Z g - c. - -* + * 1,-0 


' J=1 


s c _ 

_j rr rr 


and 


S JJ C ’ Jr + Eyr S rr c ’rr =0 


J =1 

If equation (A17) also holds, the previous equation becomes 

n 

hi c ’jr * 0 


(A19) 




It can he seen from figure 5 that if equations (A17) and (A19) 
hold, then any setting Xp will affect only Xp. Equations (A17) 

and (A18) thus are sufficient for the noninteraction conditions 
specified to hold. 
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Figure 4. - Combined general oontrolled-englne-system configuration. 



































